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Abstract. 

Recent studies of core-collapse supernovae have revealed the existence of 
two distinct classes of massive supernovae (SNe): 1) very energetic SNe (Hy- 
pernovae), whose kinetic energy (KE) exceeds 10 52 erg, about 10 times the KE 
of normal core-collapse SNe, and 2) very faint and low energy SNe (E < 0.5 

x 10 51 erg; Faint supernovae). These two classes of supernovae are likely to be 
"black-hole-forming" supernovae with rotating or non-rotating black holes. We 
compare their nucleosynthesis yields with the abundances of extremely metal- 
poor (EMP) stars to identify the Pop III (or first) supernovae. We show that 
the EMP stars, especially the C-rich type, are likely to be enriched by black- 
hole-forming supernovae. 



1. Introduction and Summary 

Stars more massive than ~ 25 M & form a black hole at the end of their evolu- 
tion. Stars with non-rotating black holes are likely to collapse "quietly" ejecting 
a small amount of heavy elements (Faint supernovae). In contrast, stars with 
rotating black holes are likely to give rise to very energetic supernovae (Hy- 
pernovae). We present distinct nucleosynthesis features of these two classes of 
"black- hole- forming" supernovae. Nucleosynthesis in Hypernovae is character- 
ized by larger abundance ratios (Zn,Co,V,Ti)/Fe and smaller (Mn,Cr)/Fe than 
normal supernovae, which can explain the observed trend of these ratios in 
extremely metal-poor stars. Nucleosynthesis in Faint supernovae is character- 
ized by a large amount of fall-back. We show that the abundance pattern of 
the recently discovered most Fe-poor star, HE0107-5240, and other extremely 
metal-poor stars are in good accord with those of black-hole-forming supernovae, 
but not pair-instability supernovae. This suggests that black-hole-forming super- 
novae made important contributions to the early Galactic (and cosmic) chemical 
evolution. Finally we discuss the nature of First (Pop III) Stars. 

2. Hypernovae and Faint Supernovae 

Type Ic Hypernovae 1998bw and 2003dh were clearly linked to the Gamma- 
Ray Bursts GRB 980425 (Galama et al. 1998) and GRB 030329 (Stanek et al. 
2003; Hjorth et al. 2003), thus establishing the connection between long GRBs 
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Figure 1. The explosion energy and the ejected Ni mass as a function of 
the main sequence mass of the progenitors for several supernovae/hypernovae 
(Nomoto et al. 2003). 



and core-collapse supernovae (SNe). SNe 1998bw and 2003dh were exceptional 
for SNe Ic: they were as luminous at peak as a SN la, indicating that they 
synthesized 0.3 - 0.5 M & of 56 Ni, and their kinetic energy (KE) were estimated 
as E 5 i = £'/10 51 erg ~ 30 (Iwamoto, Mazzali, Nomoto, et al. 1998; Woosley, 
Eastman, & Schmidt 1999; Nakamura et al. 2001a; Mazzali et al. 2003). 

Other "hypernovae" have been recognized, such as SN 1997ef (Iwamoto et 
al. 2000; Mazzali, Iwamoto, & Nomoto 2000), SN 1999as (Knop et al. 1999; 
Hatano et al. 2001), and SN 2002ap (Mazzali et al. 2002). These hypernovae 
span a wide range of properties, although they all appear to be highly energetic 
compared to normal core-collapse SNe. The mass estimates, obtained from 
fitting the optical light curves and spectra, place hypernovae at the high-mass 
end of SN progenitors. 

In contrast, SNe II 1997D and 1999br were very faint SNe with very low 
KE (Turatto et al. 1998; Hamuy 2003; Zampieri et al. 2003). In the diagram 
that shows E and the mass of 56 Ni ejected M( 56 Ni) as a function of the main- 
sequence mass M ms of the progenitor star (Figure therefore, we propose that 
SNe from stars with M ms > 20-25 M have different E and M( 56 Ni), with a 
bright, energetic "hypernova branch" at one extreme and a faint, low-energy 
SN branch at the other (Nomoto et al. 2003). For the faint SNe, the explosion 
energy was so small that most 56 Ni fell back onto the compact remnant. Thus 
the faint SN branch may become a "failed" SN branch at larger M ms . Between 
the two branches, there may be a variety of SNe (Hamuy 2003). 

This trend might be interpreted as follows. Stars with M ms < 20-25 M & 
form a neutron star, producing ~ 0.08 ± 0.03 M 56 Ni as in SNe 1993J, 19941, 
and 1987A. Stars with M ms > 20-25 M & form a black hole; whether they become 
hypernovae or faint SNe may depend on the angular momentum in the collapsing 
core, which in turn depends on the stellar winds, metallicity, magnetic fields, 
and binarity. Hypernovae might have rapidly rotating cores owing possibly to 
the spiraling-in of a companion star in a binary system. 
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Figure 2. Evolution of the central density and temperature for the 30 and 
90 M Q models (Umeda & Nomoto 2004) 



3. Synthesis of 56 Ni in ~ 100M o Stars 

The light curve modeling of the unusually bright hypernova SN1999as suggests 
that the progenitor is a core-collapse supernova and the ejected 56 Ni mass is 
as large as ~ 3 - 4M . Motivated by SN 1990as, Umeda & Nomoto (2004) 
have investigated how much 56 Ni can be synthesized in core-collapse massive 
super novae. 

The evolutions of several very massive stars with initial masses of M < 
IOOMq and low metallicity {Z = Zq/200) have been calculated from the main- 
sequence to "hypernova" explosions. The synthesized 56 Ni mass increases with 
the increasing explosion energy and the progenitor mass. Umeda &: Nomoto 
(2004) found that for the explosion energy of 3xl0 52 ergs, for example, the 56 Ni 
mass of up to 2.2, 2.3, 5.0, and 6.6 M@ can be produced for the progenitors 
with masses of 30, 50, 80 and 100 Mq, that are sufficiently large to explain SN 



Figure [2 shows the evolution of the central density and temperature for the 
30 and 9OM0 models. More massive stars have larger specific entropy at the 
center, thus having higher temperature for the same density. For 9OM , the 
evolutinary track is very close to (but outside of) the "e~e + pair-instabillity 
region" of T < 4/3 where T denotes the adiabatic index. The evolution of 
the central temperature and density is significantly different between the 30 
and 9OM models during Si-burning at T 9 = T/10 9 K =2.5 - 4. The central 
temperature and density of the 9OM model oscillate several times. This is 
because in such massive stars radiation pressure is so dominant that T is close 
to 4/3, and thus the inner core of the stars easily expands with the nuclear energy 
released by Si-burning. Once it expands, the temperature drops suddenly, the 
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central Si-burning stops, and the stellar core turns into shrink. Since only small 
amount of Si is burnt for each cycle, this pulsations occur many times. 

Umeda & Nomoto (2004) found from the study of 8O-1OOM stars that 
the number of the oscillations depends on the convective parameter larger 
fk increases the number of the oscillation. This is because for a larger fj. fresh 
Si is mixed more efficiently into the center, which increase the lifetime of this 
stage. The situation is similar to the breathing pulse phase at the end of He- 
burning, while the variation of the temperature and density is much larger in Si 
burning than in the breathing pulse phase. The amplitude of the temperature 
and density variation is larger for more massive stars, which suggests more and 
more drastic oscillations would occur for larger mass stars. 



4. Nucleosynthesis in Hypernovae 

In core-collapse supernovae/hypernovae, stellar material undergoes shock heat- 
ing and subsequent explosive nucleosynthesis. Iron-peak elements are produced 
in two distinct regions, which are characterized by the peak temperature, T pea k, 
of the shocked material. For T pea ^ > 5 x 10 9 K, material undergoes complete 
Si burning whose products include Co, Zn, V, and some Cr after radioactive 
decays. For 4 x 10 9 K < T pca k < 5 x 10 9 K, incomplete Si burning takes place 
and its after decay products include Cr and Mn (Nakamura et al. 1999). 

4.1. Supernovae vs. Hypernovae 

The right panel of Figure |3] shows the composition in the ejecta of a 25 Mq 
hypernova model (E§\ = 10). The nucleosynthesis in a normal 25 Mq SN model 
(.E51 = 1) is also shown for comparison in the left panel of Figure (Umeda & 
Nomoto 2002). 

We note the following characteristics of nucleosynthesis with very large 
explosion energies (Nakamura et al. 2001b; Nomoto et al. 2001; Umeda & 
Nomoto 2005): 

(1) Both complete and incomplete Si-burning regions shift outward in mass 
compared with normal supernovae, so that the mass ratio between the complete 
and incomplete Si-burning regions becomes larger. As a result, higher energy 
explosions tend to produce larger [(Zn, Co, V)/Fe] and smaller [(Mn, Cr)/Fe], 
which can explain the trend observed in very metal-poor stars (Umeda & Nomoto 
2005). 

(2) In the complete Si-burning region of hypernovae, elements produced by 
a-rich freezeout are enhanced. Hence, elements synthesized through capturing 
of a-particles, such as 44 Ti, 48 Cr, and 64 Ge (decaying into 44 Ca, 48 Ti, and 64 Zn, 
respectively) are more abundant. 

(3) Oxygen burning takes place in more extended regions for the larger KE. 
Then more O, C, Al are burned to produce a larger amount of burning products 
such as Si, S, and Ar. Therefore, hypernova nucleosynthesis is characterized by 
large abundance ratios of [Si,S/0], which can explain the abundance feature of 
M82 (Umeda et al. 2002). 




Figure 3. Abundance distribution against the enclosed mass M r after the 
explosion of Pop III 25 M Q stars with _E 51 = 1 (left) and E 5 i = 10 (right) 
(Umeda & Nomoto 2002). 



4.2. Hypernovae and Zn, Co, Mn, Cr 

Hypernova nucleosynthesis may have made an important contribution to Galac- 
tic chemical evolution. In the early galactic epoch when the galaxy was not yet 
chemically well-mixed, [Fe/H] may well be determined by mostly a single SN 
event (Audouze & Silk 1995). The formation of metal-poor stars is supposed 
to be driven by a supernova shock, so that [Fe/H] is determined by the ejected 
Fe mass and the amount of circumstellar hydrogen swept-up by the shock wave 
(Ryan, Norris, & Beers 1996). Then, hypernovae with larger E are likely to in- 
duce the formation of stars with smaller [Fe/H], because the mass of interstellar 
hydrogen swept up by a hypernova is roughly proportional to E (Ryan et al. 
1996; Shigeyama & Tsujimoto 1998) and the ratio of the ejected iron mass to E 
is smaller for hypernovae than for normal supernovae. 

In the observed abundances of halo stars, there are significant differences 
between the abundance patterns in the iron-peak elements below and above 
[Fe/H]~ -2.5 - -3. 

(1) For [Fe/H]< —2.5, the mean values of [Cr/Fe] and [Mn/Fe] decrease 
toward smaller metallicity, while [Co/Fe] increases (McWilliam et al. 1995; 
Ryan et al. 1996). 

(2) [Zn/Fe]~ for [Fe/H] ~ -3 to (Sneden, Gratton, & Crocker 1991), 
while at [Fe/H] < —3.3, [Zn/Fe] increases toward smaller metallicity (Cayrel et 
al. 2004). 

The larger [(Zn, Co) /Fe] and smaller [(Mn, Cr) /Fe] in the supernova ejecta 
can be realized if the mass ratio between the complete Si burning region and the 
incomplete Si burning region is larger, or equivalently if deep material from the 
complete Si-burning region is ejected by mixing or aspherical effects. This can 
be realized if (1) the mass cut between the ejecta and the compact remnant is 
located at smaller M r (Nakamura et al. 1999), (2) E is larger to move the outer 
edge of the complete Si burning region to larger M r (Nakamura et al. 2001b), 
or (3) asphericity in the explosion is larger. 
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Figure 4. Observed abundance ratios of [Zn, Co, Cr, Mn/Fe] vs [Fe/H] 
(Cayrel et al. 2004) compared with (15 M Q , E S1 = 1) and (25 M , E 51 =3Q) 
models (Tominaga et al. 2005). 



Among these possibilities, a large explosion energy E enhances a-rich freeze- 
out, which results in an increase of the local mass fractions of Zn and Co, while Cr 
and Mn are not enhanced (Umeda & Nomoto 2002, 2005). Models with E$\ = 1 
do not produce sufficiently large [Zn/Fe]. To be compatible with the observa- 
tions of [Zn/Fe] ~ 0.5, the explosion energy must be much larger, i.e., £"51 > 20 
for M > 2OM , i.e., hyper nova-like explosions of massive stars (M > 25M & ) 
with E§i > 10 are responsible for the production of Zn. 

In the hypernova models, the overproduction of Ni, as found in the simple 
"deep" mass-cut model, can be avoided (Umeda h Nomoto 2005). Therefore, if 
hypernovae made significant contributions to the early Galactic chemical evolu- 
tion, it could explain the large Zn and Co abundances and the small Mn and Cr 
abundances observed in very metal-poor stars (Fig. ^ Tominaga et al. 2005). 



5. Extremely Metal-Poor (EMP) Stars 

Recently the most Fe deficient and C-rich low mass star, HE0107-5240, was dis- 
covered (Christlieb et al. 2002, 2004). This star has [Fe/H] = -5.3 but its mass 
is as low as 0.8 Mq. This would challenge the recent theoretical arguments that 
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Figure 5. (left) The post-explosion abundance distributions for the 25 M© 
model with the explosion energy E§i = 0.3. (right) Elemental abundances of 
the C-rich, most Fe deficient star HE0107-5240 (Christlieb et al. 2004: filled 
circles), compared with a theoretical supernova yield (Umeda & Nomoto 2003, 
2005). 



the formation of low mass stars, which should survive until today, is suppressed 
below [Fe/H] = -4 (Schneider et al. 2002). 

The important clue to this problem is the observed abundance pattern of 
this star. This star is characterized by a very large ratios of [C/Fe] = 4.0 and 
[N/Fe] = 2.3, while the abundances of elements heavier than Mg are as low as 
Fe (Christlieb et al. 2002). Interestingly, this is not the only extremely metal 
poor (EMP) stars that have the large C/Fe and N/Fe ratios, but several other 
such stars have been discovered (Aoki et al. 2002). Therefore the reasonable 
explanation of the abundance pattern should explain other EMP stars as well. 
We show that the abundance pattern of C-rich EMP stars can be reasonably ex- 
plained by the nucleosynthesis of 20 - 130 M Q supernovae with various explosion 
energies and the degree of mixing and fallback of the ejecta. 

5.1. The Most Fe-Poor Star HE0107-5240 

We consider a model that C-rich EMP stars are produced in the ejecta of (al- 
most) metal-free supernova mixed with extremely metal-poor interstellar mat- 
ter. We use Pop III pre-supernova progenitor models, simulate the supernova 
explosion and calculate detailed nucleosynthesis (Umeda & Nomoto 2003). 

In Figure |S] (right) we show that the elemental abundances of one of our 
models are in good agreement with HE0107-5240, where the progenitor mass is 
25 M Q and the explosion energy £"51 = 0.3 (Umeda & Nomoto 2003). 

In this model, explosive nucleosynthesis takes place behind the shock wave 
that is generated at M r = 1.8 M and propagates outward. The resultant abun- 
dance distribution is seen in Figure [3 (left), where M r denotes the Lagrangian 
mass coordinate measured from the center of the pre-supernova model. The 
processed material is assumed to mix uniformly in the region from M r = 1.8 
Mq and 6.0 M@. Such a large scale mixing was found to take place in SN1987A 
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and various explosion models (Hachisu et al. 1990). Almost all materials below 
M r = 6.0 Mq fall back to the central remnant and only a small fraction (/ = 2x 
10~ 5 ) is ejected from this region. The ejected Fe mass is 8 x 10 -6 Mq. 

The CNO elements in the ejecta were produced by pre-collapse He shell 



burning in the He-layer, which contains 0.2 Mq 12 C. Mixing of H into the He 
shell-burning region produces 4 x 10 -4 Mq 1 N. On the other hand, only a 
small amount of heavier elements (Mg, Ca, and Fe-peak elements) are ejected 
and their abundance ratios are the average in the mixed region. The sub-solar 
ratios of [Ti/Fe] = —0.4 and [Ni/Fe] = —0.4 are the results of the relatively 
small explosion energy (£51 = 0.3). With this "mixing and fallback", the large 
C/Fe and C/Mg ratios observed in HE0107-5240 are well reproduced. 

In this model, N/Fe appears to be underproduced. However, N can be pro- 
duced inside the EMP stars through the C-N cycle, and brought up to the surface 
during the first dredge up stage while becoming a red-giant star (Boothroyd &: 
Sackmann 1999). 
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Figure 6. Averaged elemental abundances of stars with [Fe/H] = —3.7 
(Cayrel et al. 2004) compared with the hypernova yield (left: 25 Mq, E51 — 
20), and the normal SN yield (right: 15 M©, E 51 = 1). 



5.2. EMP Stars from VLT Observations 

The "mixing and fall back" process can reproduce the abundance pattern of the 
typical EMP stars without enhancement of C and N. Figure H3 shows that the 
averaged abundances of [Fe/H] = —3.7 stars in Cayrel et al. (2003) can be fitted 
well with the hypernova model of 25 Mq and £51 = 20 (left) but not with the 
normal SN model of 15 Mq and £51 = 1 (right) (Tominaga et al. 2005). 



6. First Stars 

It is of vital importance to identify the first generation stars in the Universe, 
i.e., totally metal-free, Pop III stars. The impact of the formation of Pop III 
stars on the evolution of the Universe depends on their typical masses. 

6.1. High Mass vs. Low Mass 

Recent numerical models have shown that, the first stars are as massive as ~ 100 
Mq (Abel et al. 2002). The formation of long-lived low mass Pop III stars may 
be inefficient because of slow cooling of metal free gas cloud, which is consistent 
with the failure of attempts to find Pop III stars. 
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Table 1. The results of the stability analysis for Pop III and Pop I stars. 
O and x represent that the star is stable and unstable, respectively. The 
e- folding time for the fundamental mode is shown after x in units of 10 4 yr 
(Nomoto et al. 2003). 
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If the most Fe deficient star, HE0107-5240, is a Pop III low mass star that 
has gained its metal from a companion star or interstellar matter (Yoshii 1981), 
would it mean that the above theoretical arguments are incorrect and that such 
low mass Pop III stars have not been discovered only because of the difficulty 
in the observations? 

Based on the results in the earlier section, we propose that the first gen- 
eration supernovae were the explosion of ~ 20-130 Mq stars and some of them 
produced C-rich, Fe-poor ejecta. Then the low mass star with even [Fe/H] < — 5 
can form from the gases of mixture of such a supernova ejecta and the (almost) 
metal-free interstellar matter, because the gases can be efficiently cooled by 
enhanced C and O ([C/H] ~ -1). 

6.2. Pair Instability SNe vs. Core Collapse SNe 

In contrast to the core-collapse supernovae of 20-130 Mq stars, the observed 
abundance patterns cannot be explained by the explosions of more massive, 130 

- 300 Mq stars. These stars undergo pair-instability supernovae (PISNe) and 
are disrupted completely (e.g., Umeda & Nomoto 2002; Heger & Woosley 2002), 
which cannot be consistent with the large C/Fe observed in HE0107-5240 and 
other C-rich EMP stars. The abundance ratios of iron-peak elements ([Zn/Fe] 
< -0.8 and [Co/Fe] < -0.2) in the PISN ejecta (Fig.® Umeda & Nomoto 2002; 
Heger & Woosley 2002) cannot explain the large Zn/Fe and Co/Fe in the typical 
EMP stars (McWilliam et al. 1995; Norris et al. 2001; Cayrel et al. 2003) and 
CS22949-037 either. Therefore the supernova progenitors that are responsible 
for the formation of EMP stars are most likely in the range of M ~ 20 — 130 Mq, 
but not more massive than 130 Mq. This upper limit depends on the stability 
of massive stars as will be discussed below. 

6.3. Stability and Mass Loss of Massive Pop III Stars 

To determine the upper limit mass of the Zero Age Main Sequence (ZAMS), 
Nomoto et al. (2003) analyzed a linear non-adiabatic stability of massive (80M Q 

- 300M Q ) Pop III stars using a radial pulsation code. Because CNO elements are 
absent during the early stage of their evolution, the CNO cycle does not operate 
and the star contracts until temperature rises sufficiently high for the 3a reaction 
to produce 12 C. We calculate that these stars have A"cno ~ 1-6 — 4.0 x 10~ 10 , 
and the central temperature T c ~ 1.4 x 10 s K on their ZAMS. We also examine 
the models of Pop I stars for comparison. 
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Tabled shows the results for our analysis. The critical mass of ZAMS Pop 
III star is 128Mq while that of Pop I star is 94M . This difference comes from 
very compact structures (with high T c ) of Pop III stars. 

Stars more massive than the critical mass will undergo pulsation and mass 
loss. We note that the e- folding time of instability is much longer for Pop III 
stars than Pop I stars with the same mass, and thus the mass loss rate is much 
lower. These results are consistent with Ibrahim et al. (1981) and Baraffe et al. 
(2001). However, the absence of the indication of PISNe in EMP stars might 
imply that these massive stars above 13OM0 undergo significant mass loss, thus 
evolving into Fe core-collapse rather than PISNe. 
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Figure 7. Model light curve (thick line: Nomoto et al. 2005) compared with 
the observation of SN 2002ic (filled circles; Deng et al. 2004). 

SNe 1997cy and 1999E were initially classified as Type Iln because they 
showed Ha emission. SN 2002ic would also have been so classified, had it 
not been discovered at an early epoch. SN 1997cy {z = 0.063) is among the 
most luminous SNe discovered so far {My < —20.1 about maximum light), and 
SN 1999E is also bright (M v < -19.5). Both SNe 1997cy and 1999E have been 
suspected to be spatially and temporally related to a GRB (Germany et al. 
2000; Rigon et al. 2003). However, both the classification and the associations 
with a GRB must now be seen as highly questionable in view of the fact that 
their replica, SN 2002ic, appears to have been a genuine SN la at an earlier 
phase (Hamuy et al. 2003; Deng et al. 2004). 

Nomoto et al. (2005) calculated the interaction between the expanding 
ejecta and CSM (circumstellar matter). For the supernova ejecta, the carbon 
deflagration model W7 was used. After day ~ 350, the light curve starts declin- 
ing. To reproduce the declining part of the light curve, we add the outer CSM 
of 0.2 M where the density declines sharply as n = 6. This implies that the 
total mass of CSM is ~ 1.3 Mm. 
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